Anderson Localization of Expanding Bose-Einstein Condensates in Random Potentials 
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We show that the expansion of an initially confined interacting ID Bose-Einstein condensate can exhibit 
Anderson localization in a weak random potential with correlation length <j r . For speckle potentials the Fourier 
transform of the correlation function vanishes for momenta k > 2/cr R so that the Lyapunov exponent vanishes 
in the Born approximation for k > l/<x R . Then, for the initial healing length of the condensate £i„ > tr R the 
localization is exponential, and for £i„ < a R it changes to algebraic. 
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Disorder in quantum systems can have dramatic effects, 
such as strong Anderson localization (AL) of non-interacting 
particles in random media (H. The main paradigm of AL is 
that the suppression of transport is due to a destructive inter- 
ference of particles (waves) which multiply scatter from the 
modulations of a random potential. AL is thus expected to 
occur when interferences play a central role in the multiple 
scattering process 0J. In three dimensions, this requires the 
particle wavelength to be larger than the scattering mean free 
path, I, as pointed out by Ioffe and Regel yfl. One then finds a 
mobility edge at momentum k m = 1/1, below which AL can 
appear. In one and two dimensions, all single-particle quan- 
tum states are predicted to be localized 14J,|5l|6(], although for 
certain types of disorder one has an effective mobility edge in 
the Born approximation (see Ref. J3l and below). A crossover 
to the regime of AL has been observed in low dimensional 
conductors and recently, evidences of AL have been 

obtained for light waves in bulk powders and in 2D dis- 
ordered photonic lattices ll 111 . The subtle question is whether 
and how the interaction between particles can cause dereal- 
ization and transport, and there is a long-standing discussion 
of this issue for the case of electrons in solids 11211 . 

Ultracold atomic gases can shed new light on these prob- 
lems owing to an unprecedented control of interactions, a per- 
fect isolation from a thermal bath, and t he p ossibilities of 
designing controlled random 1 13, 14, 15, Kj 17] or quasi- 



random 11811 potentials. Of particular interest are the studies of 
localization in Bose gases 11191. 120I1 and the interplay between 
interactions and disorder in Bose and Fermi gases 
Localization of expanding Bose-Einstein condensates (BEC) 
in random potentials has been reported in Refs. (HII1II1. 
However, this effect is not related to AL, but rather to the 
fragmentation of the core of the BEC, and to single re- 
flections from large modulations of the random potential in 
the tails Il5tl . Numerical calculations (l5, 23, 24] confirm 



this scenario for parameters relevant to the experiments of 
Refs. QlIIlEzl]. 



In this Letter, we show that the expansion of a ID interact- 
ing BEC can exhibit AL in a random potential without large 
or wide modulations. Here, in contrast to the situation in 



Refs. E [H, 11711 . the BEC is not significantly affected by 
a single reflection. For this weak disorder regime we have 
identified the following localization scenario on the basis of 
numerical calculations and the toy model described below. 

At short times, the disorder does not play a significant 
role, atom-atom interactions drive the expansion of the BEC 
and determine the long-time momentum distribution, T>(k). 
According to the scaling theory l25ll . T>{k) has a high- 
momentum cut-off at l/£i n , where = h/y/4mjl and \i are 
the initial healing length and chemical potential of the BEC, 
and m is the atom mass. When the density is significantly de- 
creased, the expansion is governed by the scattering of almost 
non-interacting waves from the random potential. Each wave 
with momentum k undergoes AL on a momentum-dependent 
length L(k) and the BEC density profile will be determined 
by the superposition of localized waves. For speckle poten- 
tials the Fourier transform of the correlation function vanishes 
for k > 2/ct r , where cr R is the correlation length of the disor- 
der, and the Born approach yields an effective mobility edge at 
l/er R . Then, if the high-momentum cut-off is provided by the 
momentum distribution T>(k) (for £j n > er R ), the BEC is expo- 
nentially localized, whereas if the cut-off is provided by the 
correlation function of the disorder (for £i„ < er R ) the localiza- 
tion is algebraic. These findings pave the way to observe AL 
in experiments similar to those of Refs. 115Ula . ll7ll . 



We consider a ID Bose gas with repulsive short-range in- 
teractions, characterized by the ID coupling constant g and 
trapped in a harmonic potential Va (z) = mw 2 z 2 /2. The 
finite size of the trapped sample provides a low-momentum 
cut-off for the phase fluctuations, and for weak interactions 
(n ^> mg/h 2 , where n is the ID density), the gas forms a true 
BEC at low temperatures l2rj|l . 

We treat the BEC wave function i/j(z, t) using the Gross- 
Pitaevskii equation (GPE). In the presence of a superimposed 
random potential V(z), this equation reads: 



iHdtip = 



-h 2 

2m 



d 2 + V ho (z) + V(z)+g\iP\ 2 -/i 



4>, (i) 



where tp is normalized by J'dzl'i/'l 2 = N, with N being the 
number of atoms. It can be assumed without loss of gener- 



2 



ality that the average of V(z) over the disorder, (V), van- 
ishes, while the correlation function C(z) = (V(z')V(z' +z)) 
can be written as C(z) = V^c(z/a R ), where the reduced 
correlation function c(u) has unity height and width. So, 
Vr = v (V 2 ) is the standard deviation, and a R is the corre- 
lation length of the disorder. 

The properties of the correlation function depend on the 
model of disorder. Although most of our discussion is gen- 
eral, we mainly refer to a ID speckle random potential M27I1 
similar to the ones used in experiments with cold atoms 
EH Emm. It is a random potential with a truncated 
negative exponential single-point distribution Jz" 



V[V(z)} 



exp[-(V(z) 



Vr 



Vr)/Vr] (Vjz) 



(2) 



where is the Heaviside step function, and with a correlation 
function which can be controlled almost at will II 1 711 . For a 
speckle potential produced by diffraction through a ID square 
aperture U^H^l, we have 



C(z) = V£c{z/a R ); c(u) = sin 2 (u)/u 2 



(3) 



Thus the Fourier transform of C(z) has a finite support: 

d{k)=Vl<T R c(ka R ); c(k) = y/n/2(l-n/2)e{l-K/2), (4) 

so that C(k) = for k > 2/er R . This is actually a general 
property of speckle potentials, related to the way they are pro- 
duced using finite-size diffusive plates M27I1 . 

We now consider the expansion of the BEC, using the fol- 
lowing toy model. Initially, the BEC is assumed to be at 
equilibrium in the trapping potential Vh (z) and in the ab- 
sence of disorder. In the Thomas-Fermi regime (TF) where 
jU 3> ftio, the initial BEC density is an inverted parabola, 
n(z) = Qi/g)(l - z 2 /L 2 P )Q(l - \z\/L TF ), with L TF = 



being the TF half-length. The expansion is in- 



duced by abruptly switching off the confining trap at time 
t = 0, still in the absence of disorder. Assuming that the 
condition of weak interactions is preserved during the expan- 
sion, we work within the framework of the GPE ([TJ. Repulsive 
atom-atom interactions drive the short-time (t < 1 juS) expan- 
sion, while at longer times (t 1/w) the interactions are not 
important and the exp ansion becomes free. According to the 
scaling approach [25], the expanding BEC acquires a dynam- 
ical phase and the density profile is rescaled, remaining an 
inverted parabola: 

ip(z,t) = (d[z/b{t),Q]/y/Ht)) cxp {\mz 2 b{t)/2hb{t)} 1 

(5) 

where the scaling parameter b(t) = 1 for t = 0, and b{t) ~ 
V2ujt fort > 1/u) (l5|. 

We assume that the random potential is abruptly switched 
on at a time to ^> 1/uj. Since the atom-atom interactions are 
no longer important, the BEC represents a superposition of 
almost independent plane waves: 



dk 

7^ 



ip(k, t) exp(ifcz). 



(6) 



The momentum distribution T>(k) follows from Eq. ©. For 
t ^> 1 jio, it is stationary and has a high-momentum cut-off at 
the inverse healing length l/£i n : 



2?(ft) = MM)| 5 



(1 - fc^)9(l - fc6n), (7) 



with the normalization condition J_ 00 dkT>(k) = N. 

According to the Anderson theory fc-waves will expo- 
nentially localize as a result of multiple scattering from the 
random potential. Thus, components exp(ikz) in Eq. © will 
become localized functions <f>k(z). At large distances, 4>k{z) 
decays exponentially, so that ln|</>fc(z)| ~ — 7(fc)|z|, with 
7(fc) = 1/L(k) the Lyapunov exponent, and L(k) the local- 
ization length. The AL of the BEC occurs when the indepen- 
dent fc-waves have localized. Assuming that the phases of the 
functions 4>k{z), which are determined by the local properties 
of the random potential and by the time to, are random, un- 
corrected functions for different momenta, the BEC density 
is given by 

/>oo 

n (z) ee (|V(z)| 2 ) = 2 / AkV(k){\^ k {z)\ 2 ), (8) 
Jo 

where we have taken into account that T>(k) = T>(—k) and 
(|<Mz)| 2 ) = <|0- fc (z)| 2 ). 

We now briefly outline the properties of the functions (f>k (z) 
from the theory of localization of single particles. For a weak 
random potential, using the phase formalism 12811 the state 
with momentum k is written in the form: 

4> k {z) = r{z) sin [00)] ; d z <t> k = kr(z) cos [9(z)} , (9) 

and the Lyapunov exponent is obtained from the relation 
■y(k) = — lim| z |_, 00 (log [r{z)\ /\z\). If the disorder is suffi- 
ciently weak, then the phase is approximately kz and solving 
the Schrodinger equation up to first order in \d z 9(z)/k — 1|, 
one finds 12811 . 



7 (fc) ~ (V2Tr/8a R )(VR/E) z (ka R yc(2ka R ), (10) 

where E = h 2 k 2 /2m. Such a perturbative (Born) approxima- 
tion assumes the inequality 



V R a R « (h 2 k/m)(ka R )^ 2 , 



(ID 



or equivalently -f(k) <C k. Typically, Eq. ( fTTT i means that the 
random potential does not comprise large or wide peaks. 

Deviations from a pure exponential decay of <j>k are ob- 
tained using diagrammatic methods l29ll . and one has 



7T 2 7 (fc) 



du u sinh(7ru) x 



(12) 



1 + cosh(7ra) 



exp{-2(l + W 2 ) 7 (fc)|z|}, 



where 7 (fc) is given by Eq. ( TTOb . Note that at large dis- 
tances (7(fc)|z| > 1), Eq. ([T2]i reduces to (\4> k (z)\ 2 ) ~ 

ny 2 /6Ay / 2^(k)\z\ 3 / 2 ) exp{-2 7 (fc)|z|}. 
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The localization effect is closely related to the properties of 
the correlation function of the disorder. For the ID speckle 
potential the correlation function C(k) has a high-momentum 
cut-off 2 / <r R , and from Eqs. (0]i and ( fT0T > we find 



7(fc)=7o(ft)(l-fcr R )e(l-fcoiO; 7o (*:) = 



irm 2 V^<7 R 
2h 4 k 2 ■ 



(13) 



Thus, one has j(k) > only for ka R < 1 so that there is 
a mobility edge at l/er R in the Born approximation. Strictly 
speaking, on the basis of this approach one cannot say that the 
Lyapunov exponent is exactly zero for k > l/a R . However, 
direct numerical calculations of the Lyapunov exponent show 
that for k > 1 ja R it is at least two orders of magnitude smaller 
than 7o(l/<7 R ) representing a characteristic value of -f(k) for 
k approac hing l/er R . For cr R > 1/j/m, achievable for speckle 
potentials 11711 and for Vr satisfying Eq. (fTTT i with k ~ l/cr R , 
the localization length at A: > l/er R exceeds 10cm which is 
much larger than the system size in the studies of quantum 
gases. Therefore, k = l/er R corresponds to an effective mo- 
bility edge in the present context. We stress that it is a general 
feature of optical speckle potentials, owing to the finite sup- 
port of the Fourier transform of their correlation function. 

We then use Eqs. (Q, dT2b and ( fT3l for calculating the den- 
sity profile of the localized BEC from Eq. ([8]). Since the high- 
momentum cut-off of T>(k) is 1/ £i n , and for the speckle poten- 
tial the cut-off of 7 (A;) is 1/ cr R , the upper bound of integration 
in Eq. (|8]> is /s c = min{l/£in, l/c R }. As the density pro- 
file no(z) is a sum of functions (\<t>k{z)\ 2 ) which decay expo- 
nentially with a rate 2-f(k), the long-tail behavior of Uq{z) is 
mainly determined by the components with the smallest 7(fc), 
i.e. those with k close to k c , and integrating in Eq. ([8]) we limit 
ourselves to leading order terms in Taylor series for T)(k) and 
7 (A;) at k close to k c . 

For > cr R , the high-momentum cut-off k c in Eq. ([8]) is 
set by the momentum distribution T>(k) and is equal to l/£i„. 
In this case all functions (\<fik(z)\ 2 ) have a finite Lyapunov 
exponent, 7(fc) > 7(l/£in), and the whole BEC wave function 
is exponentially localized. For the long-tail behavior of no{z), 
from Eqs. (01, (O and ( TTZb we obtain: 



n (z) cx |zr 7 / 2 cxp{-27(l/ein)|2:|}; & 



(14) 



Equation ( TBI assumes the inequality 7(l/£in)M ^ L or 
equivalently 7o(fc c )(l - cr R /£ in )|z| > 1. 

For £j n < <T R , fc c is provided by the Lyapunov expo- 
nents of (| (j>f. (z) | 2 ) so that they do not have a finite lower 
bound. Then the localization of the BEC becomes alge- 
braic and it is only partial. The part of the BEC wave func- 
tion, corresponding to the waves with momenta in the range 
l/er R < k < l/£in, continues to expand. Under the condition 
7o(fec)(l — £in/c R )M ^ 1 f° r me asymptotic density distri- 
bution of localized particles, Eqs. dS) and (fT2l yield: 



n (z) cx \z\ 



(15) 



Far tails of no (z) will be always described by the asymp- 
totic relations ( TBi i or (Q3), unless £j n = er R . In the special case 



°f £in = cr R , or for £j n very close to er R and at distances where 
7o(fc c )|(l - i 2 Jo- 2 )z\ < 1, still assuming that 7o(fc c )|z| > 1 
we find n n (z) cx |z|~ 3 . 

Since the typical momentum of the expanding BEC is 
l/£i n , according to Eq. ( fTTI ). our approach is valid for Vr <C 
^(^in/CTR) 1 / 2 . For a speckle potential, the typical momen- 
tum of the waves which become localized is l/er R and for 
£in < °r the restriction is stronger: Vr <C ^ten/cR) 2 - These 
conditions were not fulfilled, neither in the experiments of 
Refs. Jll[i3,[i7tl, nor in the numerics of Refs. |fl5j IBl I2I . 



We now present numerical results for the expansion of a 
ID interacting BEC in a speckle potential, performed on the 
basis of Eq. (Q~|i. The BEC is initially at equilibrium in the 
combined random plus harmonic potential, and the expansion 
of the BEC is induced by switching off abruptly the confining 
potential at time t = as in Refs. |H1 [Ho E3l ■ The dif- 
ferences from the model discussed above are that the random 
potential is already present for the initial stationary conden- 
sate and that the interactions are maintained during the whole 
expansion. This, however, does not significantly change the 
physical picture. 

The properties of the initially trapped BEC have been dis- 
cussed in Ref. 12211 for an arbitrary ratio £i n /<7R. For £j n <C c R , 
the BEC follows the modulations of the random potential, 
while for > cr R the effect of the random potential can 
be significantly smoothed. In both cases, the weak random 
potential only slightly modifies the density profile fl22ll . At 
the same time, the expansion of the BEC is strongly sup- 
pressed compared to the non-disordered case. This is seen 
from the time evolution of the rms size of the BEC, Az = 



2 , in the inset of Fig. [TJ At large times, the BEC 



density reaches an almost stationary profile. The numerically 
obtained density profile in Fig. [TJ shows an excellent agree- 
ment with a fit of tiq(z) from Eqs. 0, (HJ and (Tl2l . where a 
multiplying constant was the only fitting parameter. Note that 
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Figure 1: (color online) Density profile of the localized BEC in a 
speckle potential at t — 150/lj. Shown are the numerical data (black 
points), the fit of the result from Eqs. I0, <[8l» and d 1 2b [red solid 
line], and the fit of the asymptotic formula d 1 41 > [blue dotted line]. 
Inset: Time evolution of the rms size of the BEC. The parameters are 
Vr = O.lfi, Cm = O.OILtf, and a R = 0.78£ in . 



4 




&n=0-01L TF 
V R =0.05|i 



0.000 



pean Union (FINAQS consortium and grants IST-2001-38863 
and MRTN-CT-2003-505032), the ESF program QUDEDIS, 
and the Dutch Foundation FOM. LPTMS is a mixed research 
unit 8626 of CNRS and University Paris-Sud. 



1.4 1.6 1.8 2.0 



Figure 2: (color online) a) Lyapunov exponent 7 e ff in units of 1/Ltf 
for the localized BEC in a speckle potential, in the regime £i„ > 
(T R . The solid line is 7(l/£m) from Eq. d 1 3b . b) Exponent of the 
power-law decay of the localized BEC in the regime £i„ < er R . The 
parameters are indicated in the figure. 



Eq. ([8]) overestimates the density in the center of the localized 
BEC, where the contribution of waves with very small k is 
important. This is because Eq. ( fT"3l ) overestimates -f(k) in this 
momentum range, where the criterion ( fTTT i is not satisfied. 

We have also studied the long-tail asymptotic behavior of 
the numerical data. For £j n > <j r , we have performed fits of 
|z| -7 / 2 e -27crf l z l to the data. The obtained 7 e g- are in excel- 
lent agreement with 7(1/ £; n ) following from the prediction of 
Eq. (TPfl) . as shown in Fig. |2^. For £j n < er R , we have fitted 
| z |-Afr to (jjg data. The results are plotted in Fig.|2j5 and show 
that the long-tail behavior of the BEC density is compatible 
with a power-law decay with /3 e g ~ 2, in agreement with the 
prediction of Eq. < fT3T ). 

In summary, we have shown that in weak disorder the ex- 
pansion of an initially confined interacting ID BEC can ex- 
hibit Anderson localization. Importantly, the high-momentum 
cut-off of the Fourier transform of the correlation function for 
ID speckle potentials can change localization from exponen- 
tial to algebraic. Our results draw prospects for the obser- 
vation of Anderson localization of matter waves in experi- 
ments similar to those of Refs. |f5, 2 17]. For Vr = 0.2/i, 
£in = 3cr R /2 and a R = 0.27/im, we find the localization length 
L(l/£j n ) ~ 460^111. These parameters are in the range of ac- 
cessibility of current experiments fnll . In addition, the local- 
ized density profile can be imaged directly, which allows one 
to distinguish between exponential and algebraic localization. 
Finally, we would like to raise an interesting problem for fu- 
ture studies. The expanding and then localized BEC is an ex- 
cited Bose-condensed state as it has been made by switching 
off the confining trap. Therefore, the remaining small interac- 
tion between atoms should cause the depletion of the BEC and 
the relaxation to a new equilibrium state. The question is how 
the relaxation process occurs and to which extent it modifies 
the localized state. 
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